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Abstract 

We study the cosmological effect of the simple scalar-torsion (0 + ) mode in Poincare gauge theory 
of gravity. We find that for the non-constant (affine) curvature case, the early evolution of the 
torsion density px has a radiation-like asymptotic behavior of a -4 with a representing the scale 
factor, along with the stable point of the torsion pressure (Pt) and density ratio Pt/pt 1/3 in 
the high redshift regime > 0), which is different from the previous result in the literature. We 
use the Laurent expansion to resolve the solution. We also illustrate our result by the execution 
of numerical computations. 
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I. INTRODUCTION 



The recent cosmological observations, such as those from type la supernovae 



microwave background radiation 



BE 



4| , large scale structure (5 



6[ and weak lensing [7| , reveal 




cosmic 



that our universe is subject to a period of acceleration. In general, there are two ways to 
explain the phenomenon of the late-time accelerating universe 8] either by modifying the 
left- and right-handed sides of Einstein equation, called modified gravity and modified matter 
theories, respectively. For modified gravity theories, the acceleration is accounted as a part 
of the gravitational effect, while modified matter theories are constructed by including some 
negative pressure matter that could result in the expanding effect. In this study, we adopt 
the viewpoint of an alternative gravity theory by selecting so-called Poincare gauge th eory 



12 



13|- 



(PGT) |9h11|]. which is also suitable to describe the late-time accelerating behavior 

PGT is under the consideration of gauging Poincare group P4 = R 1 ' 3 x 0(1, 3) for gravity, 
and it sets out from a Riemann-Cartan spacetime (M, g,V), where M is a differentiable 
manifold, g is a metric on M, and V is a general affine metric-compatible connection with 
Vg = so that it has a canonical decomposition into V = V + K with V the Riemannian 
part and K the contortion tensor written in terms of the torsion tensor T of V. As a result, 
PGT is in general a theory of gravitation {9), [l(3] with torsion, which couples to the spin 
source. The theory comprises a degenerate case of the Einstein's general relativity (GR) of 
the vanishing torsion T and the Einstein-Cartan theory with a torsion field equation 
(TFE) algebraically coupled to the intrinsic spin of the source, resulting in a non-dynamical 
torsion field. 

In this work, we concentrate on a specific quadratic theory with only the scalar-torsion 
mode in PGT, which possesses the dynamical torsion field. This particular mode is the 
simple + mode, which is one of the six modes: 0^, 1 ± and 2^ labeled by spin and parity, 



based on the linearized theory 15 



161 ] . We remark that the mode interacts with intrinsic 



spins of fermions [12l . Il7j . However, its contribution is considered to have largely diminished 
from the early universe to the present time so that its effect in the current stage must be 
slight. On the other hand, since the + mode has no interaction with any fundamental 



source 



12 



li 



one could imagine that it remains to have a considerable portion within the 



current universe. Consequently, this mode naturally becomes the subject to study 



19]. 



In view of the FLRW cosmology, the vanishing spin current of the scalar-torsion mode 



renders a set of nonlinear equations which address the evolutions of the metric and torsion 



field. Under the positivity energy argument, Shie, Nester and Yo (SNY) in 12] have observed 
two separate cases: one has a constant affine curvature (R = const) which violates the 
positive energy condition; and the other subject to the condition forms a system of nonlinear 
ordinary differential equations (ODEs). In the former, it has a late time de Sitter space 



asymptote 



yet the torsion energy density could be negative, and the equation of state 



(EoS) of the torsion field has an interesting behavior [131 ] . In the latter, no obvious analytic 
solution is found so that numerical methods are generally applied. In particular, one of 
interesting features in our previous study [3] for the latter case is that the torsion EoS 
has an asymptotic behavior in the high redshift regime, whereas some other studies in 



the literature 



12 



20] point out that the affine curvature R, torsion scalar $ and Hubble 



parameter H = a/a are oscillatory during the cosmological evolution. Since the oscillating 



behaviors do not appear in our work [131 ]. a thorough study is clearly needed. In this paper, 
we find a proof from a semi-analytical solution in the large curvature regime to support 
our non-oscillatory result. In order to demonstrate the conformity with our semi- analytical 
solution, we will also present the numerical analysis. 

II. SCALAR-TORSION COSMOLOGY 
A. Formulation 

In this note, we explore a specific scalar-torsion mode in PGT called the simple + mode, 
given by [3] 

W = y R + y A R2 + J ( T 'i* ^ + 2T ^ ^ ~ 4Tk Tk ) > C 1 ) 
with the positive coefficients of do i and b are required by the positivity energy argument. Un- 



der the FLRW metric, the field equations have been shown in Eqs. (2.11)-(2.13) in Ref. 13 ]. 
The Friedmann equations for the scalar-torsion mode are given by 



H 



. Pc_ _ Pai + Pt 
3cto 3cto 



jj _ Pc+ Ptot _ Pc+ Pm + Pt (2) 
2a 2a U 

where ao = (87rG) _1 in GR, the subscript M represents the ordinary matter including both 
dust (m) and radiation (r), and px and Pt correspond to the torsion density and pressure 



of the effective geometric effect other than GR, defined as 



Pt = 3pH 



R 



6/i 



is v;~ ' b 

p T 



(3H - $) 5 



24 



-R 2 



PT = -(R-R) + 



respectively. From Q, we can discuss the torsion EoS, wt, defined by [13] 

Pt 



wt 



Pt 



(3) 



(4) 



B. Semi-analytical solution in high redshift 

In this context, we provide the semi-analytical solution of the positive energy scalar- 
torsion mode in the large scalar affme curvature limit i? S> Qp/b which is commonly achieved 
in the high redshift regime (a 1). In such situation, we write the energy density of ordinary 
matter and torsion in series expansion of a(t) as 

o {0) o {0) 

ym yr 

Pm = — jf H t-i 

a 6 ar 



Pt 
o (0) 

ym 



J2 A -k ak 



(5) 



k=—c 



respectively 

Before the analysis, we first follow the rescaling of the parameters in 

5o = aQ/m 2 b, Si = ai/m 2 b, i = t ■ m, p = a + a\, 
H 2 = H 2 /m 2 , $ = $/m, R = R/m 2 , 



13|, 



(6) 



where m 2 = p$ /3ao- The equations of motion [l3| can be rewritten as dimensionless 
equations: 



respectively, where \ = pi°^ / pm . In (jUJ), we have taken the approximation of R ^> 6p/b 
for the high redshift regime. With the above rescaling, we shall argue that the lowest order 
of p T does not exceed a -4 in the following discussion. We formulate the statement as a 
theorem. 

Theorem 1. In the high redshift regime (a <C 1), pr = 0(a~ 4 ). 
Proof. First we expand 

oo 

H 2 (t)=J2 r ka k ~\ (r k <oc) (11) 

k=—c 

where c is some integer, so that we have 

k=—c x ' 

Using O, ©, (HU) and flT2]), we obtain 

£=^(f>r fe a*-A+^|, (13) 



jJL J P a 



, = _3~ gl (EL C ^-4ha^)-f 

2 ai(Er=- c ^^ fe - 4 )+5 ' 

Substituting (fTTj) . (fl2|) . (fl3|) and (fT4|) into ([8]), and comparing the lowest power (requiring 

c > — 1, otherwise losing its leading position) of a in the high redshift, a <C 1, we derive the 

following relation 



2 



„2 „3 



ai 

— I c~ ■ r" 



c 2 + (5-^)c + 4 



0, (15) 

which leads to r_ c = 0ifc>las0< 5o/A < 1 and c 2 + (5 — ^)c + 4 7^ 0. This is equivalent 
to say that (ITT]) has the form 

# = ^; + ^ + ^H ^ r 4 H (16) 

a 4 a 13 ar a 

Finally, we achieve our claim from fl2]) that 

PT fX 1 V £V 2 



Pm 



X , 1 "\ , Z^o , ri r 2 r 3 \ / 1 



a 4 a 3 / V a 4 a 3 a 2 a / V a 4 



Note that the last equality follows since ro 7^ X-, which will be explained later. This is the 
end of the proof. □ 



We now write the expansion in ([IT]) , by the theorem above, simply as 

oo 

A -k ak ■ 



p 



PT 
(0) 



(18) 



We shall only take first few dominating terms for a sufficient demonstration. By the proce- 
dure in the proof of the theorem, we can as well compare terms of various orders to yield 
the following relations, 

2 



0(a~ 10 ) 
0(0 
0(a- 8 ) 



P 



0. 



2 1 



a i 



A, 



a i 
P 



P 
A- 



1 



IP 
L P 



A s A 



ai 
P 

P 



-± (A A + X ) A 2 
P 



0. 



(19) 
(20) 



a i 



- 3A 2 + ^ (A 2 (2 + 5A 3 ) 
P 



ISA, (A 4 + x)) 



0. 



(21) 



From (Tl9l . (1201) and (ED), one concludes a relation, 



A 3 = -Ji 



(«o + QiJ 
ai 



<-l, 



(22) 



with A\,A 2 and left as arbitrary constants to be determined by initial conditions and 
( |T0|) . Note that ( 122]) implies ri = —d /di < in ( fill . However, due to the observational 
data that a = 1 at the current stage, the radiation density is much smaller than the dust 
density (% <C 1), whereas the torsion density is the same order as the dust density, as seen 
from (EC 



(o) 
Pt 



[(r - x) + r 2 + •••]- (1 + N) ~ 0(1) . 



(23) 



Subsequently, we have that [(r — x) + r 2 + •••] < max{0(l), 0(|rx|)}, along with the 
assumption < 00 for each A;. As a result, we conclude that r^, for all k ^ 1, should not 
be too large, which forbids the possibility r = x- This argument shows the validity of the 
last equality in (TTTj) with the non- vanishing 0(1/ a 4 ) coefficient. 
From (J4]), via the continuity equation [13j, we obtain 



-1 



Pt 



1 / 4A 4 a -4 + 3A 3 a~ 



1 



-4, 



3p T '3 V. A A a~ A + A 3 a" 3 y 3 V" 4 4 

where the prime "/" stands for d/dliaa and we have used ([TBI for a -C 1. 



(24) 
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Wt (a) [tvt-1/31 (b) 




FIG. 1. Evolutions of (a) wt and (b) \wt — 1/3| as function of the redshift z and the scale parameter 
a, respectively, where the parameters and initial conditions are chosen as 5o = 2, 01 = 1, Hq = 2, 
R = 14 and x = pf ] / pm = 3.1 x 10~ 4 . 

C. Numerical computations 

In this subsection, we perform numerical computations to support the analysis above. As 
an illustration, we take the parameters &o = 2 and a± — 1 and initial conditions 

H( z = 0)=H = 2, R(z = 0)=R = 14, 

and show the evolutions of wt, and R in Figs. [U andEb, respectively. 



In Fig. QJi, we demonstrate the EoS of torsion as a function of the redshift z. As seen from 
the figure, in the high redshift regime wt approaches 1/3, which indeed shows an asymptotic 
behavior. Fig. [TJd indicates that \wt — 1/3| approximates a straight line in the scale factor 
a in the log-scaled coordinate since the slope in the log-scaled coordinates is nearly 1. The 
singularity in the interval [0.1, 1] corresponds to the crossing 1/3 of wt- Thus, the numerical 
results concur with our semi-analytical approximation in ( 124)) . In Fig. |2j we observe that 
the behaviors of R and $ oc 1/a 2 in the high redshift regime are consistent with the results 
in ( IT3|) and (TT^j) . given by 

R^^lA 2 a-\ (25) 
$ ~ 3H oc a" 2 , (26) 

7 



(a) $ (b) 




FIG. 2. Evolutions of (a) the rescaled affine curvature R and (b) the torsion $ as functions of the 
scale parameter a in the log scale with the parameters and initial conditions taken to be the same 
as Fig. [TJ 

respectively, where H 2 ~ (x + A 4 ) a~ A from (1TT]) . Note that from (I25j) . the behavior of 
the affine curvature i? is highly different from that of the Riemannian scalar curvature 
R = —T/ao = Pm/ao, which is proportional to 1/a 3 in both matter (dust) and radiation 
dominated eras. 



III. CONCLUSIONS 

We have investigated the asymptotic evolution behaviors of the scalar-torsion mode in 
PGT. The EoS of the torsion density has an early time stable point wt(z ^> 0) — > 1/3. This 
behavior can be estimated through the semi-analytical solution via the Laurent expansion 
in the scale factor a(t) for the torsion density p?. We have shown that there indeed exists 
the lowest degree of pt in its expansion by a -4 , corresponding to the radiation-like behavior 
in the high redshift regime. By the comparison of the next leading-order term of a~ 3 in the 
field equations, we have extracted the coefficient A3 = —fi/ai, which results in the vanishing 
of the a~ 3 term in the affine curvature R, such that R is only proportional to a -2 , consistent 
with the numerical demonstration. 
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